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摘要 : 以 控制 和 影响 褐飞虱 Nilaparvata lugens (SAD 迁 飞 过 程 的 生理 生态 参数 为 依据 ， 应 用 时 空 
分 辨 率 较 高 的 中 尺度 数值 预报 模式 一 MM4 和 三 维 轨迹 计算 方法 对 我 国盛 夏 褐 飞 乔 的 北 迁 过 程 进 
行 了 动态 模拟 。 模 拟 结果 及 其 与 实测 虫 情 资料 的 对 比 研究 表明 : (1) 我 国盛 夏 褐飞虱 北 迁 的 虫 源 
地 主要 在 22.5°N~27°N; 110°E~116°E 之 间 。(2) 空中 迁 飞 路 径 有 三 条 : 主 径 取 IVAW, Bil 
径 分 别 取 10" 和 75 方位 角 。(43) 降落 虫 汇 区 有 三 个 : EMR ARIS PRE, BIR Al ASB 
西北 、 川 东北 稻 区 和 浙 东 南 、 冰 北 稻 区 。(4) 理论 模拟 与 实际 虫 情 普 查分 析 比 较 吻 合 ， 说 明 该 模 
型 可 作为 迁 飞 害虫 灾变 机 制 研 究 的 一 个 重要 工具 。 
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4 KA Nilaparvata lugens (SAD 是 危害 水 稻 的 重要 害虫 之 一 ， 具 群 聚 性 、 远 距离 迁 飞 
性 和 为 害 突 发 性 ， 它 广泛 分 布 于 中 国 、 日 本 、 朝 鲜 半 岛 、 东 南亚 、 南 亚 、 太 平 洋 岛屿 和 澳 大 
利 亚 等 国家 和 地 区 。 在 我 国 ， 神 飞 剧 每 年 冬 春季 仅 局 限于 两 广 南部 、 台 琼 诸 岛 和 云南 南部 的 
热带 地 区 存活 与 为 害 ， 春 夏 之 交 开 始 北 迁 并 陆续 降落 与 为 害 江 南大 部 ， 盛 夏 北 迁 至 江淮 流域 
及 其 以 北 地 区 ， 降 落 并 为 害 全 国 各 水 稻 主要 生长 区 ， 夏 末 秋 初 开 始 回 迁 ， 并 继续 降落 与 为 害 
上 述 部 分 地 区 。 每 年 这 种 周期 性 的 北 迁 南 返 与 为 害 给 我 国 水 稻 生 产 造 成 了 较 大 的 危害 ， 尤 其 
是 大 发 生年 ， 稻 米 损 失 更 为 惨重 。 盛 夏 ， 褐 飞 虱 降落 与 为 害 的 区 域 最 广 、 水 稻 受 害 也 最 严重 ， 
因此 ， 研 究 我 国盛 夏 褐飞虱 北 迁 过 程 就 显得 特别 重要 。 

昆虫 三 维 迁 飞 轨迹 的 研究 始 于 1983 年 ， 当 时 美国 学 者 Demino etal. (1983) 用 Reap 
(1972) 研制 的 六 层 大 尺度 大 气 扩散 三 维 轨迹 模型 SETL1 首 次 推算 出 了 小 地 老虎 迁 入 美国 依 
Ba EM a) = Ra A), IG, Wolf el al. (1987〉 也 用 该 模式 模拟 出 了 芋 带 绿 夜 蛾 迁 入 美国 
和 加 拿 大 的 轨迹 :31。 进 入 90 年 代 以 来 ， 国 外 对 三 维 轨迹 的 研究 鲜 见 报道 。 国 内 ， 胡 继 超 等 
(1997) 首次 应 用 时 空 分 辨 率 较 高 的 中 斥 度 气象 数值 预报 模式 MM2 (MM4 的 前 身 ) 对 我 国 东 
部 秋季 回迁 的 褐飞虱 进行 了 三 维 轨迹 研究 ， 并 对 以 南京 空中 某 层 为 起 、 终 点 的 昆虫 12 h 二 、 
三 维 顺 、 逆 推 迁 飞 轨迹 作 了 对 比分 析 ， 最 后 用 雷达 观测 资料 进行 了 验证 ， 证 实 了 三 维 轨迹 比 
二 维 更 精确 、 反 映 的 迁 飞 特征 也 更 全 面 [4]。 
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前 几 年 ， 我 们 进行 了 双 雷 达 昆 虫 迁 飞 观测 、 氢 气球 空中 网 捕 采 样 、 琢 斯 高 塔 吸虫 、 田 间 
忠 情 普查 分 析 和 二 维 迁 飞 轨迹 模拟 等 研究 ， 获 得 了 许多 较 有 价值 的 结果 [5~8j。 为 了 进一步 了 
解 褐飞虱 的 迁 飞 和 灾变 规律 、 揭 示 其 大 发 生 的 生理 生态 机 制 ， 本 研究 在 总 结 前 人 工作 的 基础 
上 应 用 80 年 代 末 从 美国 引进 的 中 尺度 气象 动力 数值 预报 模式 一 MM4 和 ， 模 拟 出 了 1991 年 
7 月 发 生 在 我 国 的 两 个 典型 的 神 飞 乔 北 迁 过 程 中 昆虫 在 不 同时 次 从 不 同 地 区 迁 出 、 在 不 同 层 
次 飞行 、 到 达 不 同 降落 区 的 顺 、 逆 推 三 维 迁 飞 轨迹 144 条 ， 并 对 这 些 轨迹 的 特征 进行 了 综合 
分 析 ， 确 定 出 理论 上 褐飞虱 北 迁 的 虫 源 区 、 迁 飞 路 径 和 降落 区 ， 最 后 将 这 些 结果 与 实测 虫 情 
资料 进行 对 比分 析 ， 验 证 了 该 模型 的 精确 度 。 


1 资料 与 方法 


1.1 资料 来 源 

禄 飞 虱 虫 情 资 料 为 农业 部 全 国 病虫害 综合 治理 开放 实验 室 提 供 的 1971—1995 年 期 间 神 飞 
孔 虫 情 普 查 资料 ， 其 中 重点 分 析 了 1991 年 7 ASS 49 个 植保 站 的 神 飞 起 灯 下 诱捕 量 资 料 和 
田间 上 虫 情调 查 资料 。 

气象 资料 为 南京 气象 学 院 实习 气象 台 和 江苏 省 气象 台 提 供 的 1991 年 7 月 11 一 16 A. 
18 一 23 日 每 天 00GMT 和 12GMT 两 个 时 次 、10"N 一 S0"N 5 50E ~140°E 范围 内 600 多 个 气 
象 台 站 的 高 空 和 地 面 资 料 ， 这 些 资料 经 格式 转换 、 检 误 及 客观 分 析 后 有 270 多 个 台 站 的 气象 
资料 在 模式 所 选 网 格 区 内 真正 有 效 。 
1.2 Az 
1.2.1 MM4 模 式 简介 : 中 尺度 模式 MIM4 (Mesoscale Model Version 4.0) 是 80 年代 末 由 美国 
宾夕法尼亚 州立 大 学 PSU) 和 美国 国家 大 气 研究 中 心 (NCAR) 联合 研制 的 一 个 适合 于 有 
限 区 域 的 、 多 功能 的 、 高 分 辩 率 的 中 尺度 动力 气象 数值 预报 模式 ， 目 前 已 广泛 应 用 于 西方 国 
家 的 气象 、 环 保 等 部 门 0~14J。 在 我 国 该 模式 已 开始 普遍 应 用 于 气象 系统 的 中 期 和 短期 天 气 
预报 、 中 尺度 暴雨 试验 、 各 种 气象 要 素 和 环境 物理 量 预 报 、 区 域 气候 数值 模拟 等 !2-~141。 模 
式 中 包容 海陆 分 布 、 地 形 、 植 被 、 海 温 、 气 象 要 素 等 多 方面 的 信息 ， 兼 顾 大 、 中 、 小 不 同 尺 
度 的 天 气 系统 和 各 种 各 样 的 环境 背景 ， 对 边界 层 物理 过 程 有 较 好 的 处 理 方 案 ， 可 作 天 气 预 报 、 
传播 性 病菌 及 迁 飞 性 害虫 的 三 维 动态 模拟 、 风 载 污 染 物 的 轨迹 计算 、 降 水 分 析 、 和 气象 要 素 场 
分 析 等 ， 是 目前 世界 上 比较 先进 的 中 尺度 数值 预报 工具 之 一 ， 其 系统 运行 流程 参见 图 1。 
1.2.2 褐飞虱 三 维 迁 飞 轨迹 计算 方法 介绍 : 轨迹 分 析 既 可 用 于 污染 气象 学 中 计算 大 气 污染 物 
和 放射 性 物质 的 长 距离 输送 ， 用 于 天 和气 学 中 气 块 运行 轨迹 计算 ， 也 可 用 于 卫生 防疫 学 和 植物 
保护 学 中 计算 病菌 和 迁 飞 性 害 忠 的 传播 与 迁移 1。 本 模型 在 对 实测 气象 资料 进行 格式 转换 、 
检 误 、 客 观 分 析 、 垂 直 内 插 和 初始 化 后 ， 产 生 边 界 条 件 、 运 行 主 模式 、 诊 断 分 析 、 和 输出 了 每 
小 时 一 次 的 气象 要 素 预 报 场 。 然 后 ， 以 适宜 于 褐飞虱 迁 飞 的 生理 生态 条 件 为 前 提 ， 采 用 动力 
气象 学 中 常用 的 三 维 空 气质 点 轨迹 计算 方法 对 褐飞虱 夏季 北 迁 的 顺 、 逆 推 轨 迹 进 行 了 研究 。 
计算 步骤 是 : (1) 根据 MM4 主 模式 运行 的 结果 及 神 飞 乔 生 理 上 适应 的 迁 飞 高 度 对 地 面 到 
3 600 m 高 度 之 间 每 300 m 一 层 的 各 层次 用 拉 格 朗 日 插值 法 作 温度 、 位 势 高 度 、 经 纬 问 水 平 风 
分 量 和 垂直 速度 的 内 搬 。(2) 将 ec 坐 标 系 中 产生 的 垂直 速度 o 场 转变 为 直角 坐标 系 中 以 m/s 
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图 1 MM4 模式 运行 流程 
Fig. 1 The technological process of operating MM4 
为 单位 的 w (垂直 速度 ) 场 ， 并 将 以 气压 p (hpa) 表示 的 垂直 高 度 转换 为 以 m 为 单位 的 通用 
高 度 。(3) 从 轨迹 起 点 开始 ， 以 1 h 为 时 步 ， 连 续 24 h、32 h 或 其 它 历 时 (根据 轨迹 运行 状 
况 而 定 ) 计算 各 时 步 经 向 、 纬 向 及 垂直 方向 的 速度 和 位 移 。(4) 将 各 时 次 轨迹 点 用 折线 连接 
起 来 ， 即 得 所 求 的 轨迹 线 。 所 有 轨迹 计算 完成 后 再 将 轨迹 图 形 〈 包 括 三 维 轨迹 曲线 及 其 在 各 
平面 上 的 投影 ) 用 NCARGKS 绘图 软件 输出 打印 。 
1.2.3 福 飞 乔迁 飞 轨迹 参数 的 选择 与 拟定 : 顺 推 轨迹 计算 时 ， 轨 迹 起 点 选择 为 : (1) 以 广西 
全 州 为 实测 虫 源 地 代表 ; (2) 以 强 上 升 气流 中 心 为 理论 虫 源 地 代表 ; G) 以 上 升 气流 区 边界 
上 的 点 代表 理论 虫 源 区 边界 。 轨 迹 起 点 高 度 选择 : (1) 源 地 15 人 C 层 高 度 代表 褐 飞 恒 成 层 迁 飞 
的 上 限 ; (2) 源 地 18C 层 高 度 代表 最 适 迁 飞 层 高 度 ; (3) 源 地 22C 层 高 度 代表 成 层 迁 飞 下 限 
高 度 。 轨 迹 计 算 的 时 间 起 点 选择 : (1) 以 16 时 代表 成 层 始 迁 时 次 ; (2) 以 18 时 代表 峰 前 迁 
出 时 次 ; (3) 以 20 时 代表 蜂 后 迁 出 时 次 ;(4) 以 22 时 代表 迁 出 的 终了 时 次 。 
逆 推 轨迹 计算 时 ,轨迹 起 点 选择 ;: (1) 以 江苏 江宁 县 的 点 (31°59'N, -118°37E) 为 实测 
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FBRPD (HALASRATIBAGE RRR MED; (2) 以 长 江 中 下 游 地 区 的 强 下 沉 气 
流 区 中 心 为 理论 主 降 区 代表 ; G) 以 下 沉 区 最 北 界 代 表 宰 飞 乔 迁 飞 可 到 达 的 最 北 位 置 ， 下 沉 
区 边界 上 其 它 点 代表 褐 飞 避 理论 降落 区 边界 。 轨 迹 起 点 高 度 选 择 : G) 降落 区 15t BERRE 
飞 层 最 终 上 限 高 度 2) 降落 区 18 仿 层 代表 最 适 迁 飞 层 最 终 高 度 ; G) 降落 区 22C0ERRIL 
KERA TRAR. WAARA: G) ARA 2 WARE KR KR RARE 
次 ; (2) 以 次 日 20 时 代表 降落 的 峰 后 时 次 ; (3》 以 18 时 代表 降落 的 峰 前 时 次 ; (4) 以 16 时 
代表 初 降 《 即 初始 迁 入 〉 时 次 。 

顺 、 首 推 轨迹 的 时 步 、 总 历时 的 选择 及 迁 飞 距离 的 确定 : 两 种 轨迹 时 间 步 长 都 取 1 ho UE 
飞 总 历时 根据 轨迹 在 垂直 气流 上 升 区 和 下 沉 区 的 位 置 选 择 为 : (1) 24h; (2) 32h; G) 起 
人 时 次 至 轨迹 高 度 为 零 时 次 或 至 为 零 前 一 时 刻 ; (4》 起 始 时 次 至 轨迹 入 海 时 次 或 至 入 海 前 一 
ETZ; (5) 40 h《〈 主 要 了 解 外 界 条 件 满 足 时 ， 从 源 区 南 界 迁 出 的 神 飞 剧 单 体 到 达 最 北 界 时 的 
轨迹 历程 )。 时 步 、 历 时 选 好 后 ， 确 定 总 历时 内 每 小 时 轨迹 的 空间 位 移 总 和 为 迁 飞 距离 《以 
km 为 单位 )。 

迁 飞 速率 、 迁 飞 角 度 及 迁 飞 高 度 的 计算 以 轨迹 计算 迭代 完 后 各 小 时 三 维 风 速 矢量 的 绝 
对 值 代表 迁 飞 速率 。 以 正 北方 位 为 零 起 点 ， 顺 时 针 方向 为 正 、 逆 时 针 方 癌 为 负 ， 计 算 各 轨迹 
段 的 方位 角 ; 以 朝 水 平 正 北 方向 为 零 起 点 ， 向 上 转 为 正 ， 向 下 转 为 负 ， 计 算 经 向 仰角 ; 以 朝 
水 平 正 东方 向 为 零 起 点 ， 同 上 转 为 正 ， 癌 下 转 为 久 ， 计 算 纬 向 仰角 《三 个 角度 合成 后 为 迁 飞 
角度 )。 迁 飞 高 度 确定 时 ， 将 各 次 迭代 后 的 终了 点 即 各 轨迹 段 的 末端 点 位 置 的 无 量 纲 垂 直 坐 标 
转换 为 以 m 为 单位 的 高 度 值 就 是 迁 飞 高 度 。 


2 结果 及 分 析 


2.1 轨迹 模拟 结果 及 分 析 

本 研究 模拟 计算 了 1991 年 7 月 两 个 典型 的 神 飞 剧 北 迁 过 程 中 昆虫 种 群 在 不 同时 次 、 从 不 
同 地 区 迁 出 、 在 不 同 层次 飞行 、 到 达 不 同 降 落 区 的 顺 、 逆 推 轨 迹 共 144 条 ， 从 理论 上 推断 分 
析 了 褐飞虱 的 迁 出 忠 源 地 、 空 中 迁 飞 路 径 和 降落 区 ， 研究 了 迁 飞 速率 、 迁 飞 层 高 度 、 迁 飞 层 
厚度 、 迁 飞 角 度 、 迁 飞 历 时 和 迁 飞 距离 等 迁 飞 轨迹 参数 的 时 空 变化 规律 《 另 文 发 表 )}， 分 析 了 
不 同 轨迹 间 的 差异 及 风 场 、 温 场 、 高 度 和 天 气 背 景 等 对 福 飞 虱 迁 飞 的 影响 :364， 并 用 实测 虫 情 
资料 对 理论 研究 结果 作 了 印证 。 

图 2 是 从 模拟 出 的 144 条 褐飞虱 三 维 北 迁 轨迹 中 随机 抽样 的 三 条 ， 其 中 图 2:a 是 1991 年 
7 月 13 日 16 时 从 实测 忠 源 地 中 心 迁 出 后 在 迁 飞 层 顶 飞行 24 h 后 降落 于 下 沉 气 流 区 东北 界 的 
轨迹 ， 图 2;b 是 7 月 13 日 20 时 从 上 升 气 流 区 中 心 迁 出 后 在 最 适 迁 飞 层 中 飞行 24 h 降落 于 下 
沉 气 流 区 中 心 的 轨迹 ， 而 图 2;c 则 是 7 月 20 日 16 时 从 虫 源 区 南 界 迁 出 后 在 迁 飞 层 底 飞行 32h 
后 降落 于 下 沉 气 流 区 中 心 的 轨迹 。 这 些 轨迹 比较 逼真 地 反映 了 福 飞 乔 种 群 的 迁 飞 动态 。 
2.1.1 SEHR: 经 过 对 144 条 轨迹 的 综合 分 析 发 现 : 理论 上 主 虫 源 区 基本 上 都 是 
—4s, B122.5°N~27°N, 110° E~116°E 之 间 、 以 东西 向 为 长 轴 ， 南 北向 为 短 轴 、 面 积 约 为 
495 km X 660 km 的 和 矩形 区 域 ， 种 群 迁 出 密集 区 中 心 点 位 于 《25*21’'N，112°25'EE)， 从 该 区 迁 
出 的 概率 为 0.9167。 次 级 虫 源 区 有 两 个 : 一 个 是 辣 西 北 与 鄂 东 南 交 界 的 地 区 ， 中 心 点 在 
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(29N, 11E) 左右 ， 该 区 范围 最 小 ， 面 积 约 10 kmx20 km， 迁 出 概率 为 0.0417; 另 一 个 
是 闽 东 南 与 粤 东 北 交 界 的 地 区 ， 中 心 点 在 (27N, 1I8°E) 附近 ， 面 积 约 20 kmx20 km, Œ 
出 概率 为 0.0416 (图 略 )。 从 两 次 过 程 各 时 次 上 升 气 流 的 最 大 影响 范围 来 考虑 ， 主 虫 源 区 的 
分 布 区 域 为 : 南 界 理论 上 可 达 21"47'N， 但 实际 上 该 处 已 位 于 海上 ， 不 能 形成 迁 出 事实 ， 故 
研究 中 取 陆 地 上 上 升 气流 超过 昆虫 主动 迁 飞 速率 的 最 南 界 (22.5"N) 作为 主 虫 源 区 南 界 ; 北 
界 有 一 例 达 27"38N (统计 时 未 考虑 ) ， 其 它 都 在 27°N 以 南 ; BAD LIOIOE, KFX 116°E 
(有 一 例 达 116"11 王 ， 统 计时 未 考虑 ) (图 3:b )。 因 此 ， 由 该 模型 模拟 出 的 虫 源 区 与 根据 实际 
虫 情 资料 推 断 分 析出 的 结果 (图 3:a ) 比较 吻合 。 
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图 2 1991 ¥7 RA KARA AA det 图 3 两 次 北 迁 过 程 的 虫 源 地 、 迁 飞 路 径 
过 程 的 部 分 轨迹 和 降落 区 分 析 与 模拟 
a. 13 日 16 时 始 迁 ; b. 13 日 20 时 始 迁 ; a. 由 实测 虫 情 资 料 推断 分 析 的 结果 ; 
c. 20 日 16 时 始 迁 b. 由 MM4 模式 模拟 的 结果 
Fig. 2 Some trajectories modeled from two typical Fig. 3 Analysis and simulations for emigrating 
events of N. lugens’ migration which occurred sources, migrating path and immigrating regions 
in July of 1991 a. Emigrated at 16:00 on 13th; of BPH population for the two events a. Inference 
b. Emigrated at 20:00 on 13th; and analysis of practical investigation; 


c. Emigrated at 16:00 on 20th b. Simulations from MM4 
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2.1.2 空中 迁 飞 路 径 的 分 析 : 由 主 虫 源 区 迁 出 的 福 飞 筷 一 般 取 三 条 路 径 : (1)》 迁 飞 主 径 为 东 
北 稍 偏 北方 向 ， 即 从 源 区 出 发 经 桂 东 北 、 赣 西 ， 过 皖南 到 达 皖 中 和 江苏 境内 ， 降 落 时 可 波及 
鄂 东 、 皖 北 、 鲁 驳 南 部 、 上 海 和 浙 北 一 线 范围 内 。 该 路 径 的 轨迹 起 点 与 终点 连 线 的 平均 方位 
角 为 32" 左 右 ， 变 化 于 18" 一 54" 之 间 ， 众 值 为 33" 左 右 ， 占 总 轨迹 数 的 70.4%: (2) 第 一 副 径 
为 偏 北方 向 ， 即 从 源 区 出 发 经 湘南 、 鄂 西南 到 达 川 东北 和 哪 西北 地 区 ， 轨 迹 平 均 方位 角 为 
9"53“， 变 化 于 -15$" 一 1$" 之 间 ， 众 值 为 10" 左 右 ， 占 总 轨迹 数 的 18.8%: (3) 第 二 副 径 为 偏 东 
方向 ， 即 从 源 区 出 发 经 更 南 、 粤 北 ， 到 达 浙 东南 和 闽 北 地 区 ， 轨 迹 平均 方位 角 为 7$"36， 变 
化 于 60°~85°21' Ziel, PTS AA, ABB 10.8% (图 3d. 
2.1.3 迁 入 降落 区 的 分 析 ; 模拟 发 现 ， 在 144 条 从 理论 虫 源 地 迁 出 、 取 不 同 路 径 完成 迁 飞 历 
程 的 神 飞 起 轨迹 中 ， 不 管 南端 的 起 点 在 何 处 ， 但 北端 的 落 点 都 分 布 在 三 个 密集 区 内 。 其 中 主 
降 区 在 长 江 中 下 游 稻 区 ， 即 20°-N~36°N, 115°E~122°E 之 间 ， 且 轨迹 北端 落 点 在 该 区 内 以 
(32°22'N, 118°37E) APDEA BIH, PEMA A 107/144; 第 一 副 降 区 为 川 东北 和 
哪 西北 稻 区 ， 即 32"N 一 34"N，110"E 一 114" 眉 之 间 ， 轨 迹 落 点 密集 中 心 在 〈32"19'N， 
112*29'E)， 降 落 在 该 区 的 概率 为 19/144; 第 二 副 降 区 为 浙 东 南 和 间 北 稻 区 ， 即 27"N 一 29"N， 
117°30'E~121°30’E 之 间 ， 落 点 区 中 心 在 〈28"16'N，120?58' 下 )， 在 该 区 降落 的 概率 为 18/ 
144 (Kl 3:b). 
2.2 理论 研究 结果 与 实测 虫 情 资 料 的 对 比分 析 

首先 ， 从 虫 源 地 研究 来 看 ， 理 论 上 主要 虫 源 区 在 22.$"N 一 27"N，110" 下 一 116 开 之 间 ， 迁 
出 起 点 密集 分 布 在 以 点 〈25?21'N，112?25 开 ) AP DON ARRAN BK RMA, KWER 
区 在 23°N~27°N, 110°E~115°E Ż 0, FAP DAZE (25°56'N, 111°07 ED) AABA— 
致 《 图 3)， 只 是 理论 源 区 稍 大 些 。 两 虫 源 中 心 纬 癌 偏 差 144.33 km, 4 64.22 km, E 
线 距离 为 137.86 km。 在 所 计算 的 144 条 轨迹 中 ， 有 129 条 南端 在 实测 主 虫 源 区 内 ， 有 7 了 7 条 
离 实测 主 源 区 外 缘 不 足 70 km， 另 外 8 条 可 能 因 散 迁 从 “次 级 忠 源 区 ” 迁 出 ， 可 见 吻合 率 达 
90% 以 上 。 

其 次 ， 从 迁 飞 路 径 和 降落 区 域 来 看 ， 模 拟 出 的 主 径 、 第 一 、 第 二 副 径 在 144 条 轨迹 中 除 
8 条 “ 散 迁 ”轨迹 外 其 余 均 与 由 多 年 虫 情 资 料 判 断 分 析 的 结果 基本 一 致 〈 图 3)， 吻 合 率 达 
94.44%。 降 落 区 域 与 忠 情 资料 分 析 基 本 吻合 ， 两 个 副 降 区 外 缘 与 实测 完全 一 致 ， 理 论 主 降 区 
(包括 波及 区 东 、 西 、 南 三 界 与 实际 相合 ， 北 界 比 实测 推断 的 稍 北 些 ， 纬 度 上 相差 43”， 即 
90 km 左右 。 实 测 主 降 区 最 大 概率 中 心 在 《31 57’N，11831’E)， 理 论 主 降 区 相应 中 心 在 
(32"22'N，118*37'E)， 两 者 直线 距离 相差 70.28 kmo 

上 述 对 比分 析 表 明 : MM4 是 一 种 比较 理想 的 轨迹 模拟 工具 ， 同 时 笔者 在 本 研究 中 对 不 同 
迁 飞 参数 的 选择 也 比较 合理 。 


3 结论 与 讨论 
本 研究 应 用 MM4 模式 ， 引 进 控制 和 影响 褐飞虱 迁 飞 过 程 的 生理 生态 参数 ， 结 合 三 维 轨 


迹 计 算 方 法 模拟 得 到 了 1991 年 7 月 发 生 在 我 国 局 部 地 区 的 两 个 典型 福 飞 剧 北 迁 过 程 的 三 维 轨 
迹 144 条 ， 分 析 研 究 后 发 现 ， 
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(1) Re KAP Bite tS 22.5°N~27°N, 110°E~116°E 之 间 约 500 km X 650 
km 的 区 域内 。 迁 出 后 空中 飞行 路 径 有 三 条 : 主 径 取 东北 偏 北方 向 ， 水 平方 位 角 32* 左 右 ; 第 
一 副 径 取 正 北 偏 东 方向 ， 方 位 角 10* 左 右 ; 第 二 副 径 取 正 东 偏 北方 向 ， 方 位 角 75* 左 右 。 褐 飞 
改 种 群 迁 飞 后 的 降落 区 域 有 三 个 ， 主 降 区 为 长 江 中 下 游 稻 区 ;第 一 副 降 区 为 鄂 西北 和 川 东 北 
稻 区 ;第 二 副 降 区 为 浙 东 南 和 国 北 稻 区 :为 害 的 重 灾区 在 30"N 一 34"N > 116°E~120°E 之 间 。 

(2) 盛夏 神 飞 乔迁 出 迁 入 的 初始 时 刻 都 在 16 时 左右 ， 高 峰 时 刻 在 18 一 20 时 之 间 ， 终 目 
时 刻 为 22 时 左右 。 

(3) MM4 模式 是 一 个 时 空 分 辩 率 较 高 的 中 尺度 动力 数值 预报 模式 ， 其 水 平 格 距 可 在 30 
到 150 km 之 间 调 放 ， 垂 直方 向 从 地 面 到 100 hpa 之 间 可 分 成 若干 层 ， 在 时 间 上 可 分 辩 30 s 到 
36 h 之 间 各 种 时 间 间 隔 的 气象 要 素 和 大 气 物理 量 的 变化 状况 。 三 维 轨迹 计算 方法 是 反映 大 气 
质点 或 空中 风 载 浮游 物 移动 状况 及 其 时 变 特征 的 有 效 方法 之 一 。 应 用 它们 作 昆 虫 迁 飞 轨迹 模 
拟 时 ， 只 要 选择 合适 的 初始 迁 飞 参数 和 运行 限制 条 件 就 可 以 获得 一 些 令 人 满意 的 研究 结果 。 
本 研究 理论 模拟 与 实测 推断 分 析 的 逼近 ， 正 好 说 明了 这 一 点 。 

(4) 本 研究 在 宰 飞 乔迁 出 与 迁 入 时 刻 、 上 升 与 降落 机 制 、 成 层 的 空间 范围 等 方面 考虑 了 
它 在 生理 上 的 主动 适应 性 ， 其 合理 性 有 待 进一步 探讨 。 














参考 文 献 (References) 


[1] Reap RM. An operational three-dimensional trajectory model. J. Appl. Met. 1972, 11: 1193~1202 

[2] DemioR Pet al. Spring weather pattern associated with suspected black cutworm introduction to Iowa. J. Environ. Ent. 
1987, 12: 1863~1871 

[3] Wolf WW, Westbrook J K, Sparks AN. Relationship between radar entomological measurements and atmospheric struc- 
ture in south Texas during March and April 1982. In: Sparks A N. ed. Long-range Migration of Moths of Agronomic Im- 
portance to the United States and Canada: Specific Examples of Occurrence and Synoptic Weather Patterns Conducive to Mi- 
gration» Washington: USDA/ARS-43. 1987, 98~104 

[4] 胡 继 超 ， 程 极 盖 和 等， 褐飞虱 在 我 国 东 部 地 区 秋季 回迁 的 三 维 轨迹 研究 ， 中 国 农业 气象 ，1997，18 (5): 1~6 

[5] AZK. 昆虫 雷达 应 用 初探 ， 江苏 省 首届 青年 学 术 年 会 论文 集 . 北京 : 中 国 科学 技术 出 版 社 ，1992: 347~352 

Lo] 程 遐 年 等 .褐飞虱 在 中 国 东部 秋季 回迁 的 雷达 观察 .南京 农业 大 学 学 报 ，1994,，17 G): 24~32 

[7] 程 极 盖 ， 包 云 季 等 ， 葛 东 稳 飞 虱 暴发 的 轨迹 分 析 ， 中 国 农业 气象 ，1994，15 (1): 2~5 

[L8] 程 极 益 ， 包 云 轩 等 .褐飞虱 迁 飞 轨迹 研究 ， 南 京 农 业 大 学 学 报 ，1995，18 (3): 60~67 

[9] AnthesR A, Hsie E-Y, Kuo Y-H. Description of the Penn ./NCAR Mesoscale Model Version4 (MM4) . NCAR TECH. 
Notes 1987, NCAR/TN-282+ STR. 1~65 

[10]  Karyampudi V M and Carlson T N . Analysis and numerical simulation of the Saharan air layer and its effect on easterly 
wave disturbances. J. Atmos. Sci. 1988, 45 (21): 3 102~3 136 

[11] Westphal D L; Toon © B and Carlson T N. A case study of mobilization and transport of Saharan dust. J. Atmos. Sci. 
1988, 45 (15): 2145~2175 

[12] 沈 A, B20. FARR KAT RAR GZ. RAFIR 1996, 32: 216~223 

[13] ŽA, = lg, 改进 的 MM4 模式 暴雨 实例 对 比试 验 . 气象 科学 ，1996 (2): 103~109 

[14] “ 刘 春 涛 ， 程 刨 生 ， 黑 风暴 的 沙 尘 形 成 与 输送 参数 化 及 中 尺度 数值 试验 . ARF, 1997, 55 (6): 726 一 739 

[15] ，” 朱 艳 秋 ， 雇 洞 捅 ， 三 维 空气 质点 轨迹 的 计算 ， 应 用 气象 学 报 ，1992，3 (3): 328 一 333 

[16] fat. Bik, WBE. 盛夏 褐飞虱 北 迁 大 发 生 的 气象 背景 : 个 例 分 析 . 南京 农业 大 学 学 报 ，1999，22 (4): 
35 一 40 





第 2 期 Baits: 中 国盛 夏 褐 飞 乔 北 迁 过 程 的 动态 数值 模拟 183 


Dynamical and numerical simulations on the processes of 
Nilarparvata lugens long-distance migration northward 
during the midsummer in China 


BAO Yun-xuan!, CHENG Ji -yi!» CHENG Xia-nian'» CHU Chang-shu?» SHEN Tong-li? 
(1. Nanjing Agricultural University, Nanjing 210095; 2. Nanjing Meteorological Institute. Nanjing 210044) 


Abstract: Based on the two migration events of Nilaparvata lugens (Stal) population, which had led to 
its serious infestation on rice in Changjiang valley, the main rice-growing regions in midsummer of 1991, 
and introduced into certain physiological and ecological parameters controlling the migrations a dynamical 
and numerical prediction model with higher spatio-temporal detectablity———MM4 and a means of calculat- 
ing for 3-D trajectory were used to model the dynamic state of their migrating northward . The simulations 
and the comparisons between the simulations and the practical pests investigation showed: (1) The emi- 
grating resource was an area from 22.5°N to 27°N on latitude and from 110°E to 116°E on longitude. (2) 
There were three spatial migratory paths, which consisted of a main path with an azimuth of about 32° and 
two assistant paths of about 10° and 75° respectively on azimuth. (3) There were three immigrating and 
descending regions with a main region situated in the middle and lower reaches of Changjiang valley, an as- 
sistant region joined northwestern Hubei with northeastern Sichuan and another assistant region enclosed 
southeastern Zhejiang and northern Fujian. (4) The conformity between the simulations and the practical 
pests investigation just illustrated that this model is an important means for the exploration of migratory 


pests’ catastrophic mechanism. 


Key words: Nilaparvate lugens (Stal) migration: dynamical and numerical simulations: 3-D trajectory 


analysis: migration parameters: catastrophic mechanism 


